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ABSTRACT 

Context. Sulphur is an important, volatile a element but its role in the Galactic chemical evolution is still uncertain and more obser¬ 
vations constraining the sulphur abundance in stellar photospheres are required. 

Aims. In this work we derive the sulphur abundances in red giant branch (RGB) stars in three Galactic halo globular clusters (GC) 
that cover a wide metallicity range (-2.3 <[Fe/H]< -1.2), namely M4 (NGC6121), M22 (NGC6656), and M30 (NGC7099). The 
halo field stars show a large scatter in the [S/Fe] ratio in this metallicity span, which is inconsistent with canonical chemical evolution 
models. To date, very few measurements of [S/Fe] exist for stars in GCs, which are good tracers of the chemical enrichment of their 
environment. However, some light and a elements show star-to-star variations within individual GCs and it is yet unclear whether the 
a element sulphur also varies between GC stars. 

Methods. We used the the infrared spectrograph CRIRES to obtain high-resolution (R ~ 50000), high signal-to-noise (SNR 200 per 
px) spectra in the region of the S I multiplet 3 at 1045 nm for 15 GC stars selected from the literature (6 stars in M 4, 6 stars in M 22 
and 3 stars in M 30). Multiplet 3 is better suited for S abundance derivation than the more commonly used lines of multiplet 1 at 
920 nm, since its lines are not blended by telluric absorption or other stellar features at low metallicity. 

Results. We used spectral synthesis to derive the [S/Fe] ratio of the stars assuming local thermodynamic equilibrium (LTE). We find 
mean [S/Fe]LTE = 0.58 ± 0.01 ± 0.20dex (statistical and systematic error) for M4, [S/Fe]LTE = 0.57 ± 0.01 ± 0.19dex for M22, 
and [S/Fe]LTE = 0.55 ± 0.02 ± 0.16dex for M30. The negative NLTE corrections are estimated to be in the order of the systematic 
uncertainties. We do not detect star-to-star variations of the S abundance in any of the observed GCs with the possible exception of 
two individual stars, one in M 22 and one in M 30, which appear to be highly enriched in S. 

Conclusions. With the tentative exception of two stars with measured high S abundances, we conclude that sulphur behaves like a 
typical a element in the studied Galactic GCs, showing enhanced abundances with respect to the solar value at metallicities below 
[Fe/H] -l.Odex without a considerable spread. 
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1. Introduction 


Sulphur is one of the less studied a-elements (produced by se¬ 
quential capturing of He nuclei) in stars. In general, the produc¬ 
tion site of the a-elements such as O, Ne, Mg, Si, S, Ar, Ca, 
and Ti is mainly associated with the eruptions of type II super¬ 
novae (SNe II). In contrast, the major fraction of the iron-peak 
elements abundance comes from SNe la, which contribute little 
or no a-elements. The different timescales of the occurrence of 
SNe la or II makes the [a/Fe] ratio a powerful tool for diagnos¬ 
ing the chemical evol ution and star formation history (SFH) of 

latteucci & Brocatoll990l: 

Way, the metal poor stars 
([Fe/H] < -l.Odex) form a plateau at [a/Fe] ~ 0.4dex. With 
the onset of SNe la, at metallicities above [Fe/H] ~ -l.Odex, 
the [g/Fe] ratio st arts to gradually drop reaching solar values 
(lMcWilliamlll997l) . In dwarf galaxies that experienced slower 
star formation, lower values of [a/Fe] are observed even at low 


lllg LllC LllCmiLill CVUl UliUll ailU blUL lUU 

any stellar population (|Tinslevll979tlMa 

iMcWilliani et alJl2013l) . In the Milky 
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metallicities, owing to the lower e nrichment rate of the inter¬ 
stellar medium from SNe II (see Shetrone_eLaL 2()0ll l2003[ 


IVenn et ^l2()04t iTolstov et al.ll2009[ iHendricks et al.l 2014 and 
references therein). Although different a elements are produced 
on similar timescales, they show element-to-element scatter due 
to different production mechanisms, either through He burning 
in the cores of massive stars, or during the SNe explosions them¬ 
selves. In particular, sulphur (along with Si, Ar, and Ca) is pro¬ 
duced during b oth the O-shell burning and explosive oxygen 
burning phases (iLimongi & Chief^l2003[) . 


Sulphur is an espec ially interesting g -element because it is 
not depleted onto dust (lUeda et al.ll2005h and is thus a genuine 
tracer of the "true" S abundance of the ISM and the stellar photo¬ 
spheres. In fact, su lphur is a preferred tracer of the a abund ance 
in interstellar gas (iGarnetlil 19891: ISavage & Sembachlll996h and 
in high r edshift star forming environments such as damped Ly-g 
system s (ICenturion et al.ll2000l : iNissen et'^l2004 iBowen et^ 
l2005h . The existing sulphur measurements in halo field stars 
show that the [S/Fe] ratio is increasing with decreasing [Fe/H] at 
higher metallicities, reaching an average value of about 0.4 dex 
but its behaviour at lower metallicities is still rather controver¬ 
sial. Different studies report either a bimodality of the [S/Fe] 
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ratio with a flat plateau and a branch of increasing abundance, 
reac hing [S/Fe] ~ 0.8 dex, or a plateau with a large scatter (see 


e.g. Nissen et al.l 12004 ICaffau et al. 

2005al; Nissen et al. 2007t 

ISnite et al.ll201 1; Llonsson et al. 2011 

; Matrozis et alJ 20131 and 


references therein). 


Globular clusters (GCs) represent the oldest stellar popula¬ 
tions in the Milky Way and are examples of very dense envi¬ 
ronments. Long time considered as simple stellar populations, 
nowadays we recognize their complex SFH through precise 
abundance analysis of a variety of chemical elements in individ¬ 
ual cluster member stars. All GCs, studied to date, present sig¬ 
nificant spreads and certain unique correlations in their light and 
a element abundances - generally taken a s a clue for multipl e 
enrichment phases (see the recent review bv iGratton et al.l2012l) . 
Still, sulphur abundances have been derived in very few GCs due 
to the difhculties that such measurements entail. Sulphur lines in 
the optical spectral range are generally very weak in giant stars. 
At the same time, these brighter stars are often the only accessi¬ 
ble tracers in remote and faint GCs. To date, the strong multiplet 
1 (920 nm) has been used to obtain sulphur abundances in GCs 


(Caffau et a 

J2005bHSbordone et al.l2009l:lKoch & CaffaiJ201 It 

ICaffau et al 

|2014|) but these lines are often blended with telluric 


absorption lines and one of them is positioned in the blue wing of 
the Paschen ( line and thus cannot always be reliably measured. 
This makes the Mult. 3 at 1045 nm particularly suited to measure 
the sulphur abundance. Even if the lines are not as strong as the 
components of Mult. 1, a big advantage is that there are no tel¬ 
luric features in this wavelength range. Observing the 1045 nm 
sulphur lines leads to reliable results even at lower metallicities 
and will give independent and possib ly more accurate sulphur 
abundances in GCs than using Mult. 1 dCaffau et al.ll2007l 120101 : 
iJonsson et al.ll20ll1) . 

The metal rich GC Ter7 ([Fe/H] = -0.32 dex) is the first 
GC, for which S abundances were ev er derived. The three 
stars measured by ICaffau et alJ (l2005bh are also the only ex- 
tragalactic stars with known [S/Fe] ratio, since Ter 7 is asso- 
ciated with t he disru pted Sagittarius dwarf spheroidal galaxy. 
ICaffau et alJ (l2005bll determined mean S abundance in Ter 7 
slightly lo wer (by ~ 0.1 dex) tha n in Galactic stars at similar 
metallicity. ISbordone et al.l (l2009l) analysed 9 stars in the GC 
47 Tuc ([Fe/H] = -0.72dex) and 4 stars in NGC 6752 ([Fe/H] = 
-1.54 dex), denying [S/Fe] = 0.18 + 0.14 and [S/Fe] = 0.49 + 
0.15 dex for both GCs, respectiyely, consistent with the results 
of field stars. Curiously. ISbordone et alJ (l2009h reported a star- 
to-star yariation in the obseryed [S/Fe] ratio in 47 Tuc, which 
although not significant with respect to the large measurement 
uncertainties, showed a strikingly clear correlation with the Na 
abundance. The authors considered different possibilities for the 
occunence of such correlation including that S may be inyolyed 
in the self-enrichment p rocesses of GCs throug h the ^*P(p, y/^^S 
proton capture reaction. iKoch & Caffaul (1201 1 ) analysed a single 
star in the metal poor GC NGC 6397 ([Fe/H] = -2.0dex), de¬ 
termining [S/Fe] = 0.52 ± 0.2 dex, consiste nt with the halo file d 
population. We also note the recent work by lCaffau et al.l (l2014ll . 
which proyides [S/Fe] measurements for seyeral metal rich open 
clusters and includes the GC M 4, analysed in the present work. 

Here, we present the homogeneous analysis of 15 bright 
RGB stars in three different GCs that coyer together a wide range 
of metallicities (-2.3 < [Fe/H] < -1.2 dex). Such strategy can 
proyide answers to seyeral basic questions. It is interesting to 
find out whether a dichotomy or a large spread of the [S/Fe] ra¬ 
tio, as obseryed in the halo field stars, could also be found in 
GCs. The generally large uncertainties of the [S/Fe] measure¬ 
ments in the halo field stars are mostly driyen by uncertainties 


of the atmospheric parameters (iMatrozis et al.ll2()T^ and thus, 
tracing the [S/Fe] ratio with GCs will lead to more accurate re¬ 
sults due to the much better constrained distances, grayities, and 
effectiye temperatures of GC stars from colour-magnitude dia¬ 
grams (CMD). Additionally, our data will help to check if there 
are star-to-star yariations of the S abundance in the different GCs 
and whether they correlate with the other a and light element 
abundances, such as O, Na, Si, Mg, and Al. 

The three selected clusters are amongst the closest to the Sun 
that coyer an optimally broad metallicity range and are all well 
studied objects. The proximity and simple CMDs of M4 and 
M30 make them desirable candidates for many obseryational 
and theoretical studies testing the current stellar eyolutionary 
models. M 30 is also notable for its yery low metallicity [Fe/H] = 
-2.3 dex. M22, on the other hand, is amongst the most massiye 
GCs. It has a yery complex CMD with prominent multiple popu- 
lations and a large metallicity spread (-1. 9 < [Fe/H] < -1.6 dex 
iDa Costa et alJlfoO^ iMarino et alJl2009l) . It is often considered 
as a remainder of a disrupted dwarf galaxy and will thus add 
another, important comparison object to our study. 

In this paper we present the largest compilation of sulphur 
abundance measurements in Galactic clusters to date. 


2. Target selection and observations 


We chose our targets from existing high-resolution suryeys of 
these GCs in the yisible band; M 4 stars were chosen from 
iMarino et al.l (l2008[i: M 22 stars from | Marino et al.l (I2()()^l201 ll) : 
M 30 stars from Cmetta et al.l (l2009allbh . Our intention was to 
select bright RGB stars that belong to different populations in 
the selected GCs. To this end, we targeted stars with different Na 
and O abundances in M 4 and M 30. M 22, on the other hand, 
is known to harbour stars with a wide metallicity spread, so we 
selected stars with different [Fe/H]. The adopted stellar param¬ 
eters of the obseryed stars from the literature are presented in 
Table [1] Choosing stars for which high-resolution spectroscopy 
was already ayailable is also important for constraining their at¬ 
mospheric parameters, such as effectiye temperature, grayity, mi¬ 
croturbulence yelocity, and metallicity, necessary for abundance 
determination. This is only possible if a large set of Fe I and Fe 
II lines are ayailable, which is not the case with the narrow waye- 
length coyerage of our spectra. Note, howeyer, that the effectiye 
temperatures and grayitie s of the three s t ars in M 30 were deter¬ 
mined photometrically by lCarretta et al.l (l2009bl) and likely haye 
larger uncertainties. We still opted to rely on t hem since these 
parameters were used by ICarretta et al.l (l2009bli to deriye abun¬ 
dances of other chemical elements in these stars, with which we 
aim to compare our S abundance estimates. 

To compare both approaches, we compared t he stellar pa¬ 
ramete rs fo r ~ 80 stars in comm on between the IMarino et al.l 
(l2008h and ICarretta et al.l (l2009bll stu dies. The effectiye tem¬ 
peratures determined by ICar retta et aP (l2009bl) are by ~ lOOK 
lower on ayerage than those in Marino et al.l ( 20081) for M 4 with 
cr — 45 K and the log g yalues used by ICarretta et ^ (l2009bl) 
are by 0.3 dex lower on ayerage with cr = 0.15 dex. Marino et al.l 
(I 2 OO 8 I) find a lower spread in microturbulence and a larger spread 
in meta llicities among the common stars, while ICarretta et al.l 
(l2009bl) find the opposite, i.e., a larger spread in microturbulence 
and a lower spread in metallicity. The ayerage differences are 
0.13 km s"' in microturbulence and 0.13 dex in [Fe/H] between 
the two studies. 

Our spectra were taken in seryice mode under Programme ID 
09l.B-0171(A) with the VLT cr yogenic high-resolu tion infrared 
echelle spectrograph tCRIRES iKaeufl et al.l l2004l) mounted at 
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Table 1. Atmospheric parameters and [S/Fe] ratios of the targeted stars. 


ID 

Tf// 

[K] 

logg 

[dex] 

mic 

km s“* 

[Fe/H] 

[dex] 

[Na/Fe] 

[dex] 

[0/Fe] 

[dex] 

[Mg/Fe] 

[dex] 

[Ca/Fe] 

[dex] 

V ^ 

^ mac 

km s ' 

[S/Fe]LTE^ 

[dex] 

19925 

4050 

1.20 

1.67 

-1.02 

M4' 

0.51±0.04 

0.28±0.04 

0.43±0.06 

0.19±0.03 

6.5±0.1 

0.57±0.02 

21191 

4270 

1.60 

1.60 

-1.06 

0.51±0.04 

0.34±0.04 

0.55±0.06 

0.24±0.03 

4.8±0.1 

0.58±0.02 

27448 

4310 

1.57 

1.58 

-1.12 

0.11 ±0.04 

0.51 ±0.04 

0.50±0.06 

0.29±0.03 

4.9±0.1 

0.57±0.01 

28103 

3860 

0.50 

1.62 

-1.08 

0.17±0.04 

0.50±0.04 

0.41±0.06 

0.16±0.03 

6.3±0.1 

0.65±0.02 

34006 

4320 

1.67 

1.61 

-1.06 

0.44±0.04 

0.25±0.04 

0.52±0.06 

0.26±0.03 

5.7±0.1 

0.50±0.02 

36215 

4300 

1.59 

1.53 

-1.11 

0.18±0.04 

0.48±0.04 

0.52±0.06 

0.26±0.03 

5.5±0.2 

0.59±0.02 

200005 

4000 

0.05 

2.02 

-1.94 

M222 

-0.02±0.03 

0.40±0.04 

0.46±0.04 

0.24±0.02 

11.9±0.2 

0.78±0.01 

200025 

4100 

0.67 

1.80 

-1.62 

0.26±0.03 

0.49±0.04 

0.44±0.04 

0.42±0.02 

8.5±0.2 

0.57±0.02 

200031 

4300 

0.77 

1.55 

-1.85 

0.31 ±0.03 

0.18±0.04 

0.34±0.04 

0.20±0.02 

9.4±0.2 

0.51±0.02 

200051 

4260 

0.90 

1.60 

-1.63 

0.73±0.03 

-0.05±0.04 

0.44±0.04 

0.40±0.02 

6.8±0.2 

0.58±0.02 

200061 

4430 

1.05 

1.70 

-1.78 

-0.17±0.03 

0.34±0.04 

0.54±0.04 

0.28±0.02 

12.1±0.6 

0.61±0.03 

200068 

4500 

1.30 

1.52 

-1.84 

-0.01±0.03 

0.45±0.04 

0.23±0.04 

0.22±0.02 

9.4±0.3 

0.56±0.02 

5783 

4463 

1.17 

2.32 

-2.33 

M30^ 

0.57±0.02 

0.16±0.12 



8.7±0.7 

0.57±0.03 

10849 

4365 

0.96 

2.34 

-2.40 

0.65±0.04 

-0.29±0.20 



11.9±0.9 

0.86±0.04 

11294 

4258 

0.41 

2.14 

-2.37 

0.07±0.07 

0.28±0.08 

0.51 ±0.20 


12.2±0.4 

0.53±0.02 


Notes. IDs, atmospheric parameters and abundances from iMarino et al] ll2008h . The adopted Solar values are as follows: A(Fe) = 7.48, 

A(Na) = 6.32, A(0) = 8.83, A(Mg) = 7.55, A(Ca) = 6.39. _ 

IDs, atmospheric parameters and abundances from iMarino et alj i20Q9h . The adopted Solar values are as follows: A(Fe) = 7.48, A(Na) = 6.31, 
A(0) = 8.83, A(Mg) = 7.54, A(Ca) = 6.39. 

IDs, atmospheric parameters and abundances from ICarretta et al.l (I2009al l5l. The adopted Solar values are as in iGratton et al.l (l2003h : 
A(Fe) = 7.54, A(Na) = 6.21, A(0) = 8.79, A(Mg) = 7.43, A(Ca) = 6.27. 

The macroturbulence velocity (v„„c) is a free parameter i n the fitting proced ure. 

The adopted Solar values are A(Fe) = 7.50, A(S) = 7.16 dCaffau et alj|201 Ih . 


the Nasmyth focus A of the UTl telescope (Antu). We used 
the 0.4" slit, which provides a resolving power R ~ 50000 
and a dispersion scale on the four parallel, 512px CCD chips 
of 0.05 A px '. The S I multiple! 3 was centred on chip 3. We 
used relatively short exposure times (between 30 and 300 s) with 
a nodding cycle between each exposure. The total integration 
time was chosen as to achieve a SNR ~ 200 per px in the re¬ 
gion of the sulphur lines. In the nodding technique, each star is 
moved along the slit and exposed at two different positions. This 
allows for proper sky subtraction and removal of systematic ef¬ 
fects, such as amplifier glow. The final, combined spectra are the 
mean of all short exposures. The science targets were also used 
as guide stars for the adaptive optics system. The observing log 
is presented in Table |3 


For the data reduction, we used the standard CRIRES 
pipeline (version 2.3.1) provided by ESO. The pipeline allows 
dark current subtraction, bad pixel correction, sky subtraction, 
flat-held correction, non-linearity correction, wavelength cali¬ 
bration from a Th-Ar lamp, and an extraction and combination 
of the final ID spectrum from the science images. We noticed 
that the wavelength calibration was suboptimal due to the very 
few Th-Ar lamp lines with sufficient brightness available in the 
observed region, so we additionally refined the wavelength so¬ 
lution using the available stellar lines in the spectra. Einally, the 
spectra were continuum normalised by fitting the continuum re¬ 
gions with a high order polynomial. 


Table 2. Observing log. 


ID 

V 

[mag] 

Date 

Exp. time 
[s] 

SNR 

[px-'] 

19925 

11.04 

M4 

Jun. 23/24 2013 

8x60 

220 

21191 

11.70 

Sep. 13/14 2013 

8x 120 

230 

27448 

11.73 

Sep. 13/14 2013 

8x 120 

230 

28103 

10.71 

Apr. 5/6 2013 

6x60 

240 

34006 

11.87 

Sep. 6/7 2013 

8x 120 

160 

36215 

11.80 

Sep. 6/7 2013 

8x 120 

140 

200005 

10.92 

M22 

Jun. 19/20 2013 

8x60 

230 

200025 

11.52 

Jun. 19/20 2013 

12x60 

190 

200031 

11.64 

Jun. 19/20 2013 

8x 120 

180 

200051 

12.04 

Jun. 20/21 2013 

10 X 120 

140 

200061 

12.21 

Aug. 9/10 2013 

14 X 180 

100 

200068 

12.30 

Aug. 10/11 2013 

16 X 120 

170 

5783^ 

12.71 

M30 

Jul. 27/28 2013 

8x300 

140 

5783^ 

12.71 

Aug. 22/23 2013 

8x300 

140 

10849 

12.55 

Jul. 27/28 2013 

12 X 300 

100 

11294 

12.09 

Aug. 22/23 2013 

6x300 

220 


Notes. This star was observed in two nights. The stated SNR is for 
the combined spectrum of all 16 x 300 s exposures. 


3. Results 

3.1. Spectral synthesis 

We derived the S abundance for each of the observed stars 
through a full spectral synt hesis using the LTE cod e “Spec¬ 
troscopy Made Easy” (SME lWlenti & Piskunoviri996ll . We de- 


Article number, page 3 of|9] 


























A&A proofs: manuscript no. niki6_v2 



X (nm) 

Fig. 1. Synthetic spectra illustrating the relative contribution of the Fe I 
line (blue dash-dotted line) at 1045.5 nm to the S I line at the same wave¬ 
length (red line) with respect to the combined spectral feature (thick 
dashed black line) at the metallicities of the three GCs assuming Solar 
abundance mix. 

rived our line list from the VALD database (iKupka et al.lll99^ 
I 2 OOOI) . The transition parameters of the lines in multiple! 3 and 
the Fe I line at 1047.0 nm are presented in Tabled There is no se¬ 
vere line blending in the wavelength region of multiple! 3. Only 
the S I line at 1045.3 nm is blended with a Fe I line but the latter 

has a marginal contribution at [Fe/H]-1 dex and practically 

disappears at lower metallicities. The relative contribution of this 
Fe I line lying on top of the S I feature is shown in Figure [T] for 
the metallicities of the three GCs. The atomic data for this blend¬ 
ing iron line is also given in Table [3] All other transitions on top 
of the S triplet are of a negligible strength but we note that we use 
the complete line-list of the VALD database, including molecu¬ 
lar bands in the respective wavelength interval in our synthesis. 
The region is free of telluric absorption lines. We interpolated 
the new grid of Kuruc43 plane-parallel, one-dimensional mod¬ 
els without convective overshoot. These incl ude the g-enhanced 
opaci ty distribution functions (AODFNEW; ICastelli & Kurucj 
l2003ln . The models fully cover the parameter space of all stars 
in our sample with steps AT,,// = 250 K, A log g = 0.5 dex, and 
A[Fe/H] = 0.5 dex. Kurucz models have also been used for the 
derivation of the stellar parameters (Table[T]). We further discuss 
the impact of different atmospheric models to the S abundances 
in Section 3.2. 

SME performs a minimisation algorithm iterating a de- 
hned set of free parameters. In our case, we used the [S/Fe] ratio 
and the macroturbulence velocity of the stars as free parameters 
to simultaneously ht the three S lines from multiple! 3 and the 
only prominent Fe line in the region (see Figures |2] [2 lH). The 
stellar parameters, as presented in Table [1] were kept constant. 
The goodness of ht of the Fe I line shows the suitability of the in¬ 
put atmospheric parameters, since it is only affected by changes 
in the broadening parameter - the macroturbulence velocity in 
this case. 

The estimated [S/Fe] ratio in LTE and the best-ht values for 
the macroturbulence velocity for each star together with the sta¬ 
tistical uncertainties are presented in the last two columns of 
Table [1] We used the value of A(S) = 7.16dex and A(Fe) = 


' http://kurucz.harvard.edu/grids.html 
^ http://wwwuser.oats.inaf.it/castelli/ 


7.50 dex as in ICaffau et al.l (1201 ih for the Solar abundance of 
sulphur and iron, respectively. 

We hnd very homogeneous [S/Fe] ratios for all 6 stars in 
M4 with a mean [S/Fe]LTE = 0.58 + 0.01 dex with an intrin¬ 
sic spread ctq = 0.04 dex. The intrinsic spread is calculated as 
ctq = ^fs/Fe]“ ^ where cr[s/Fe] is the standard devia¬ 

tion of the [S/Fe] ratio from the mean of all stars and < Srand > 
is the average random error of the individual measurements. The 
random errors (erand) are the formal uncertainties based on the 
goodness of the fit- Five of the observed stars in M 22, re¬ 
gardless of their metallicity, have consistent abundances with 
a mean [S/Fe]LTE = 0.57 + 0.01 dex with an intrinsic spread 
(To = 0.03 dex. There is only one star (ID 200005) that shows 
somewhat higher [S/Fe]LTE = 0.78 + 0.01 dex. The situation is 
similar with M 30. Two stars show consistent [S/Fe] ratios with 
a mean [S/Fe]LTE = 0.55 +0.02dex and the third one (ID 10849) 
has an abnormally high [S/Fe]LTE = 0.86 + 0.04 dex. Note, 
however, that the derived [S/Fe] ratios are strongly dependent 
on the effective temperatures and moderately on the gravities, 
and metallicities of the stars (see Section 3.2). Thus, reasonable 
changes in the effective temperatures of the stars with abnormal 
S abundances can bring them in line with the other stars in the 
sample, e.g. an increase of Tgff by 100 K for star 200005 in M 22 
and by 200 K for star 10849 in M30. Star 200005 also has a 
somewhat lower [Fe/H] estimate ~ 0. 1 dex lower than the mean 
of the metal-poor population of M 22 (iMarino et al.ll200^ . As it 
is indicated in Table 0] increasing the metallicity of the star by 
0.1 dex will bring the [S/Fe] ratio down by 0.1 dex, already in a 
reasonable agreement with the rest of the stars in M 22. 

The reference Fe I line in our spectra for star 200005 is also 
not well ht by the synthetic spectrum, with ax^ estimation about 
2.5 times higher than the average for the other stars and cannot be 
improved by only changing the macroturbulence velocity. This 
means that the model spectrum does not reproduce correctly the 
equivalent width of the line, which is another indication that the 
stellar parameters for this star might be off. We are able to obtain 
very good hts to the iron line by modifying the effective tempera¬ 
ture, gravity, and metallicity of the star but due to the degeneracy 
between these parameters, it is difficult to put any hrm conclu¬ 
sions on its sulphur abundance. The situation with star 10849 
in M 30 (it also has the lowest SNR spectrum in our sample) is 
more critical. As our comparison of the stellar par a meter d eter¬ 
mination using the methodology of ICarretta et al.l (i2009bh and 
iMarino et al.l ( 20081) for the common stars in M 4 showed, there 
are systematic differences but the standard deviations in Tgff and 
logg appear to be small. We do not have solid reasons to believe 
that the relative estimate of the effective temperature of a star in a 
homogeneous study could be off by as much as 20 0 K and it does 
not pre sent odd abundance of any other element (iCarretta et al.l 
l2009hh . 

3.2. Systematic uncertainties 

In order to estimate the systematic errors of the [S/Fe] ratio 
caused by uncertainties in the atmospheric parameters of the 
stars, we chose one star from each GC and varied its Tgff by 
±100K, logg by ±0.1 dex, Vmic by ±0.1 kms^', and [Fe/H] by 
±0.05 dex. These values are very similar to the uncertainties of 
the atmospheric parameters given in the literature. The results 
are presented in Tabled In the last column, we list the total sys¬ 
tematic error, which is calculated by combining in quadrature 
the uncertainties caused by different parameters. We note, how¬ 
ever, that this is an upper limit of the systematic uncertainties due 
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Table 3. Line parameters from VALD. 


Element 

A 

[nm] 

A" 

[eV] 

log.?/ 

Damping parameters 
Radiative Stark Waals 

Lande factor 

SI 

1045.5449' 

6.8600' 

+0.250' 

8.950' 

-5.370' 

-7.333' 

L25(P 

Fel 

1045.5404^ 

5.3930^ 

-0.905-' 

8.480^ 

-3.960-' 

-7.520^ 

1.130^ 

SI 

1045.6757' 

6.8600' 

-0.447' 

8.950' 

-5.370' 

-7.333' 

2.000' 

SI 

1045.9406' 

6.8600' 

+0.030' 

8.950' 

-5.370' 

-7.333' 

1.750' 

Fel 

1046.9653^ 

3.8840^ 

-1.187^ 

6.830'' 

-6.150^ 

-7.820^ 

1.170" 


Notes. Line data references: *^^ IZerne et ^ jl997h <^> lKuruczl j2004li <^> lBarklem etaDll200(lh <^) |Kuruczl ll2007h 


Table 4. Systematic uncertainties of [S/Fe]. 


GC Star ID ^T^ff Alogg Av^/c A[Fe/H] total 




+ 100K 

-lOOK 

+0.1 dex 

-0.1 dex 

+0.1 kms * 

-0.1 kms ' 

+0.05 dex 

-0.05 dex 


M4 

27448 

-0.18 

+0.20 

+0.04 

-0.05 

-0.01 

0.00 

-0.03 

+0.03 

±0.20 

M22 

200051 

-0.18 

+0.18 

+0.05 

-0.06 

0.00 

0.01 

-0.05 

+0.03 

±0.19 

M30 

11294 

-0.14 

+0.16 

+0.05 

-0.04 

0.00 

0.01 

-0.02 

+0.03 

±0.16 



Fig. 2. Observed (black) and synthetic (red) spectra for the stars in M 4. 
The most prominent spectral lines are marked with vertical dotted lines. 
The spectra are shifted vertically by an arbitrary value. 


to cov ariances of the different parameters (see lMcWilliam et al.l 
[1991. 

The choice of the atmosphere structure also has a marginal 
impact on the derived abundances. For instance, using a Solar- 
scaled Kurucz atmosphere model (ODFNEW) would lead to 
0.12 dex lower on average S abundances for the stars in M4. 



The differences are less pronounced in the more metal-poor 
GCs where the Solar-scaled models would render S abundances 
0.07 dex and 0.04 dex lower for M 22 and M 30, respectively. We 
also tested the Marcs 2012 spherically symmet ric models, which 
area-enhanced for metallicities below -1 dex (iGustafsson et al.l 
l2008h . These models lead to S abundances 0.07 dex higher on av¬ 
erage than the ones stated in Table[T]for all three GCs. All model 
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Fig. 5. The S multiplet 3 region in the Sun and the benchmark RGB 
star HD 83240 fitted using the same technique as for the GC stars. In 
the Sun, we estimate A(S )lte = 7.26 dex opposed to the adopted So¬ 
lar abundance of A(S) = 7.16 dex. For HD 83240 we find AfSfLXE = 
7.13 dex. Symbols are as in Figurej^ 


grids are equally densely populated in the parameter space of 
interest. 


3.3. NLTE and 3D effects 

The S lines at 1045 nm are affected by relatively large NLTE 
effects. By extrapo lating the NLTE results from the grid of 
iTakeda et al.l (l2005h to the stellar parameters of our stars, we 
found Af^LTE = -0.20dex, = -0.15dex, and Ai^ue - 

-0.10 dex for the GCs M4, M22, and M30, respectively. The 
corrections become smaller for lower metallicities. This pro¬ 
duces a weak trend of slightly increasing [S/Ee]NLTE ratios with 
decreasing metallicity with respect to the very homogeneous 
mean [S/Fe]LTE values for the three GCs stated in Section 3.1. 
We find [S/Fe]NLTE = 0.38 dex, 0.42 dex, and 0.45 dex for the 
GC M 4, M 22, and M 30, respectively - more in-line with the 
nominal [a/Fe] = 0.4 dex value for the old halo stellar popu¬ 
lation. Considering the large systematic uncertainties, the small 


increase of the [S/Fe]NLTE ratios with decreasing metallicity is 

probably insignificant. _ _ 

According to calculations b v Llonss on et al.l (1201 ih. b ased 
on the 3D model atmospheres of ICollet et al. ( 20071 200^ . the 
3D corrections for S abundances derived from multiplet 3 are 
constant and positive at all metallicities for stars with similar 
parameters to ours and a re in the order of A 3 D = - 1 - 0.20 dex. 
ICaffau et al.l (l2007l 1201 Oh explored the 3D effects of multiplet 3 
in dwarf stars and also found positive corrections based on the 
CO^BOLD 3D model atmospheres. Although it is inconsistent 
to apply their results to giant stars, there is an agreement that 
the 3D corrections for the lines at 1045 nm are always positive. 
In general, the NLTE and 3D corrections for our stars roughly 
compensate each other and are in the order of the systematic 
uncertainties. One, however, needs to be conscious when apply¬ 
ing 3D-LTE corrections to ID-NLTE corrected results since the 
abundances derived from full 3D-NLTE analysis might differ 
substantially (see Section 3.4). For this reason, and to be con¬ 
sistent with literature [S/Fe] measurements of field and cluster 
stars, we consider further in the discussion only the [S/Fe]NLTE 
results stated above. 


3 . 4 . Sulphur abundance of the Sun and HD83240 

Finally, the [S/Fe] ratio is relative to the Solar abundance of 
sulphur and therefore it is vital for the proper analysis. We de¬ 
cided to apply the same spectral sy nthesis technique to the same 
region of the Kurucz Solar atlas (iKurucz et al.lll984ll and de¬ 
rive the LTE Solar abundance of sulphur estimated from mul¬ 
tiplet 3 (Figure^. We found A(S) = 7.26 + 0.01 dex, which is 
by 0.10 dex higher than the nominal value of A(3) = 7.16 dex 
(ICaffau et al.ll201 1 I) adopted in this work. Our re sult is, how¬ 
ever, c onsistent with the Solar value estimated by ICaffau et al.l 
(l2007h using the same lines of multiplet 3 (A(3) = 7.30 dex). 
This finding will systematically decrease our derived abundances 
in LTE in the three GCs by 0.1 dex. Although the ID-NLTE and 
3D corrections for the Solar abundance derived from multiplet 
3 are small and com pensating each other - Anttf. = -0.07 d ex 
(iTakeda et alJl2005l) and A^o = -1-0.07 dex (ICaffau et al.ll2007h . it 
is likely that a NLTE analysis performed in a full 3D synthesis 
will bring the abundance down to the nominal value. The nom¬ 
inal Solar value is estimated from the S triplet at 675.7 nm and 
the [S I I line at 1082 nm, which are virtually unaffec ted by NLTE 
effects (ICaffau et al.ll2007t ICaffau & Ludwi3l2007h . 

In order to test the reliability of our results for RGB GC stars 
we also carried out an identical analysis for a well studied gi¬ 
ant benchmark sta r - HD 83240. We use d a CRIRES spectrum 
from the library of iLebzelter et al.l (1201 2lPl HD 83240 is classi¬ 
fied as a K1 III star and has a Tefr = 4680K, logg = 2.45 dex, 
Vmic = 1.3 km s ' and[Fe/H] = -0.02 dex. Alth ough, no sulphur 
abundance w as explicitly deriyed for this star, iMishenina et al.l 
(l2006il200^ measured the abundances of yarious light, Fe-peak, 
and n-capture elements and found that the star has a Solar abun¬ 
dance pattern with a slight under-abundance of carbon and oyer- 
abundance of nitrogen, consistent with a post first-dredge-up 
eyolutionary stage. Thus, we assume that HD 83240 also has So¬ 
lar abundance of sulphur. We measured [S/Fe]LTE = -0.01 dex 
(Figure|5]l. The NLTE correction for this star is expected to be in 
the order of -0.15 dex and the 3D correction is positiye and in 
the same order of magnitude. 

The Fe I line that we use to adjust the broadening in our 
spectra is well reproduced by our synthesis in both the Sun and 

^ http: //WWW. uni vie. ac. at/crirespop/index. html 


Article number, page 6 of[9] 















































































N. Kacharov et al.: Sulphur in GCs 


HD 83240. Thus we confirm the suitability of the atomic data for 
this line. 


4. Implications on the Galactic chemical evolution 
of sulphur 

We should be careful when directly comparing the stellar [S/Fe] 
abundance ratios from different studies because they are based 
on different methods, different S lines, and different assumptions 
in their final results. It has been shown that even after proper 
NLTE treatment, the various S features could re sult in offsets 
in the derived abundances (see Llonsson et alJl201 ih . Having said 
that, in Figure |6] we have compiled a large collection of known 
[S/Fe] ratios (references are provided in the figure legend and 
caption) for Galactic and cluster stars in order to roughly track 
the evolution of sulphur. All measurements have been shifted 
to the Solar value of sulphur (A(5) = 7.16dex) adopted in 
this work and NFTF corrections are applied where applicable. 
Most of the NFTF corrections come from the calculations of 
iTakeda et al. ([200^, except for the extremely metal poor stars 
of lSnite et al.l (l2Ql ill who compute their own NFTF corrections. 
We note that, although the NFTF results by both teams are com¬ 
puted for different types of stars, there are some differences in 
the overlapping regime. We did not apply 3D corrections in the 
results plotted in Figure|6]in order to be consistent with most of 
the literature discussions. 

There are S abundances derived within the entire span of 
metallicities covered by stars in Galactic clusters, although the 
sample is still limited and further studies are desirable. In all 
analysed clusters, sulphur seems to behave as a typical a el¬ 
ement, meaning that the stars in metal poor Galactic GCs 
([Fe/H] < -l.Odex) have enhanced [S/Fe] abundance ratios and 
open clusters of solar metallicity have solar [S/Fe] ratios, con¬ 
sistent with the results for field stars and chemical evolutionary 
model s. The only exceptions are Tr5 d iscussed in ICaffau et al.l 
(l2014l) and Ter 7 (ICaffau et alJ 1200^ with anomalously low 
[S/Fe] ratios. The latter is associated with the Sagittarius dwarf 
spheroidal galaxy and its low [S/Fe] ratio is consistent with the 
lower star formation efhciency of its parent galaxy. While the 
spread of sulphur on the metal poor plateau as traced by clus¬ 
ter stars is smaller on average than the spread traced by Galactic 
field stars we found one star in M 22 and one in M 30 that are 
compatible with extreme S abundances. 

Me ssier 4 is the only G C analysed by two teams indepen¬ 
dently (ICaffau et al.l 12014. and this work) and a more carefu l 
comparison between the results is possible. (ICaffau et alj|2014|) 
deriye S abundances for 10 stars using line profile fitting of mul- 
tiplet 1 at 920 nm taking also into account the telluric absorption 
lines. Although, both studies seem to agree on consistent [S/Fe] 
abundance ratios for this GC (see Figurelhli. lCaffau et alJ (120141) 
assume lower met allicity for their st ars, adopting [Fe/H] = 
-1.36 dex for M 4 ([Monaco et al.ll20 12|), w hile the mean [Fe/H] 
for our stars is = -1.08 dex dMarino et al.ll2008h . This discrep¬ 
ancy results i nto much lower ya lue of A(S) = 6.11 + 0.04 dex 
estimated by ICaffau et akl (l2014l) . compared to our estimate of 
A{S) = 6.46 + 0.02dex in NFTF. The NFTF correc tions for both 
multi plets are similar and in the order of 0.2 dex (iTakeda et al.l 
|2005|) . Both studies agree, howeyer, that there is not any signifi¬ 
cant spread in [S/Fe] amongst the RGB stars in M4. 

Finally, we explore the behayiour of sulphur with respect to 
other light and a elements. One of the goals of this study is to 
probe whether star-to-star yariations within an indiyidual GC 
are possible as is the case for elements like C, N, O, Na, Al, 


and to some extent Mg and S i, which are all corr elated between 
each other (see the reyiew by lGratton et al.ll2012L and references 
therein). In Figure Q (left panel), we plot existing measurements 
of [S/Fe] ys. [Na/Fe] for yarious GCs. In M4, M22, and M 30, 
we haye deliberately targeted stars with different Na abundances 
and Figure |7] clearly shows that there is not a S-Na correlation 
between stars from these GCs. While the high S star in M 30 is 
also Na-rich, the high S star in M22 has a low Na abundance. 
So far, wi th the only tentatiye ex ception of 47 Tuc, already dis¬ 
cussed in ISbordone et al.l (l2009h . it seems that eyen if present, 
star-to-star yariations of the sulphuric abundance in GCs are not 
caused by the same self-enrichment processes responsible for the 
formation of p-capture element abundance yariations. We note, 
howeyer, that we still lack an extensiye study that includes a large 
number of stars in a single GC. 

In the middle panel of Figure|7]we show the behayiour of the 
[S/Ca] ratio in GCs - two a elements caryed through the same 
production channels. The mean ratio is ~ 0.1 dex and there is no 
trend with metallicity, which suggests that the two elements do 
eyolye in a lock-step, as expected. There is more scatter in the 
[S/Mg] ratio, although the mean yalue is 0 dex (Figure El right 
panel). Both elements haye slightly different production chains, 
so the larger scatter is not surprising. In the pre-supernoya eyolu- 
tion. Mg is produced in the C-conyectiye shell burning, while S 
- in the O-conyectiye shell burning. During the explosiye nucle- 
osynthes is. Mg is produced in C a nd Ne burning, while S - in O 
burning (iFimongi & Chieffil 120031) . Magnesium could also par¬ 
ticipate in p-capture reactions through the MgAl hydrogen burn¬ 
ing cycle at high temperatures and hence star-to-star yariations 
of this element are prominent in some GCs. While in a majority 
of the obseryed GC stars coyering a large range of metallici¬ 
ties, the behayiour of sulphur with respect to other a elements 
fayours a typical cr-capture origin of this element, the detection 
of two yery S-rich stars (if confirmed real) may rise the question 
of additional production channels. 


5. Conclusions 

Due to the lack of suitable lines to measure sulphur in the yisi- 
ble band at low metallicities, the S abundance measurements in 
GC stars are quite scarce to date. The present work helps to sub¬ 
stantially increase the sample of cluster stars with determined S 
abundances. We measured the [S/Fe] abundance ratio from mul- 
tiplet 3 at 1045 nm in 15 RGB stars from three GCs (M 4, M 22, 
and M 30) obseryed with the high-resolution, infrared spectro¬ 
graph CRIRFS, mounted at the VFT. The multiplet 3 is suited 
for S abundance deriyation since the three lines are relatiyely 
strong eyen at low metallicity and are free of significant stel¬ 
lar blends and telluric contamination. The selected GCs coyer a 
large range of metallicities (-2.3 < [Fe/H] < -l.Odex) and the 
indiyidual RGB targets were specially selected to haye different 
Na & O abundances or eyen different Fe abundances in the case 
of the massiye GC M 22. This ensured that we haye stars belong¬ 
ing to different cluster populations and allowed us to search for 
possible star-to-star yariations of the S abundance in GCs. 

We find that the S abundances deriyed from multiplet 3 are 
yery sensitiye to changes in the effectiye temperature of the ob¬ 
seryed stars and the uncertainty of this parameter dominates the 
uncertainty of the measured abundances. The lines of multiplet 
3 are also prone to large departures from FTF and are sensi¬ 
tiye to horizontal fluctuations of the stellar atmosphere, inferring 
large 3D model corrections. The NFTF and 3D cotTections are 
similar in magnitude and roughly compensate each other, how¬ 
eyer the lack of full NFTF analysis performed within the frame- 
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★ M4 (this work) 

M4 (Caffau et al. 2014) 

★ M22 (this work) 

★ MSO (this work) 

O NGC6397(Koch & Caffau 2011) 

O NGC6752 (Sbordone et al. 2009) 

> 47 Tuc (Sbordone et ai. 2009) 

□ Ter 7 (Caffau et al. 2005) 

NGC5822 (Caffau et al. 2014) 

A NGC2477 (Caffau et al. 2014) 

+ Tr 5 (Caffau et al. 2014) 

-4 -3 -2-10 1 

[Fe/H] 

Fig- 6- NLTE corrected [S/Fe] abundance ratios as a function of m etallicity in Galactic disk and halos stars (black dots: ICaffau et alJl2005^ : 
iJonsson et al.ll201 l| : lTakeda & Takada-Hidaill2Ql IbISpite et alj[201 Ih and cluster stars (colourful symbols). All abundances are referenced to the 
adopted Solar value of A(S) = 7.16dex. The random errors of our GC stars are comparable to the symbol sizes. The maximum systematic error 
due to uncertainties of the stellar parameters is shown in the figure. 
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Fig. 7. A comparison of the S abundances with the abundances of other light and a elements in GCs: [S/Fe] vs. [Na/Fe] (left panel)’, [Fe/H] vs. 
[S/Ca] (middle panel)’, [Fe/H] vs. [S/Mg] (right panel). Different symbols indicate stars from different clusters as indicated in the legend of Figure 
1^ The maximum systematic errors are also indicated. 





work of a 3D atmospheric model introduces additional system¬ 
atic uncertainty. Nevertheless, we perceive very small star-to- 
star scatter in the majority of the observed stars in the three GCs 
consistent with a narrow plateau ([S/FejcrE = 0.57 + 0.01 dex, 
ctq = 0.03dex; [S/Fe]NLTE = 0.40dex, cro = 0.05dex), as ex¬ 
pected for the behaviour of a typical a element at low metal¬ 
licity in the Milky Way. Sulphur appears to vary in a lock step 
with Ca, a well studied a element most likely produced through 
the same nuclear channels. The mean [S/Mg] ratio is Solar but 
with a larger scatter due to the more complex production path 
of Mg, which could also participate in the self-enrichment pro¬ 
cesses that take place in GCs. 

One star in the GC M 22 and one star in the GC M 30 are 
found to be exceptionally S rich, being consistent with the high 
S measurements in some field stars (e.g. lCaffau et al.l2005ah . We 
consider the possibility that these are erroneous measurements 
due to some deviation from the true atmospheric parameters of 
the two stars but we cannot undoubtedly exclude the possibility 
of observing genuine S-rich stars in GCs. The S-rich stars do 
not show clear correlations with p-capture elements like Na and 
O, so we conclude that S is not produced during hot hydrogen 
burning reactions responsible for the formation of multiple pop¬ 
ulation in GCs. Discovering stars with confirmed extreme [S/Fe] 
ratios in different stellar environments would require the inven¬ 
tion of new production mechanisms for sulphur. 


Clearly more observations are needed, in particular large ho¬ 
mogeneous samples, to conclusively constrain the S abundances 
in GCs. 
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